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Visualizing Spacetime Curvature via Gradient Flows II: An Example of the 

Construction of a Newtonian analogue 
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This is the first in a series of papers in which the gradient flows of fundamental curvature invariants 
are used to formulate a visualization of curvature. We start with the construction of strict Newtonian 
analogues (not limits) of solutions to Einstein's equations based on the topology of the associated 
gradient flows. We do not start with any easy case. Rather, we start with the Curzon - Chazy 
solution, which, as history shows, is one of the most difficult exact solutions to Einstein's equations 
to interpret physically. We show that the entire field of the Curzon - Chazy solution, up to a region 
very "close" to the the intrinsic singularity, strictly represents that of a Newtonian ring, as has 
long been suspected. In this regard, we consider our approach very successful. As regrades the 
local structure of the singularity of the Curzon - Chazy solution within a fully general relativistic 
analysis, however, whereas we make some advances, the full structure of this singularity remains 
incompletely resolved. 

PACS numbers: 04.20.Cv, 04.20.Ha, 02.70.-c 
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I. INTRODUCTION 

Perhaps the most familiar example of an exact solution 
to Einstein's equations, generated by way of a Newtonian 
"analogue" , is the Curzon - Chazy solution wherein the 
Laplacian for a point mass in a fictitious Euclidean 3 - 
space is used to generate an exact static axially sym- 
metric vacuum solution of Einstein's equations. Unfor- 
tunately, the Curzon - Chazy solution bears little resem- 
blance to a point mass. Indeed, its singularity structure 
appears, at present, to be as complicated as that of the 
Kerr solution. 

In this paper, following the introduction given by one 
of us [1] , we consider the gradient fields of the two non- 
differential invariants of the Curzon - Chazy solution. 
These alone reveal previously unknown non-local prop- 
erties of the solution. Further, following the procedure 
given in [l| for the construction of strict Newtonian ana- 
logues, based on gradient flows, we suggest that a pure 
Newtonian ring is "almost" a complete analogue of the 
Curzon - Chazy solution. Whereas our main aim here is 
the construction of the analogue, we have had to do a 
fair amount of analysis of the Curzon - Chazy solution 
itself. Most of this analysis is relegated to Appendices. 



II. THE CURZON-CHAZY METRIC 

The line element for a static and axially symmetric 
gacetime can be written in Weyl's canonical coordinates 

ds 2 = ~e 2U dt 2 + e- 2U (e 2 ^ (dp 2 + dz 2 ) + p 2 d<f> 2 ) , (1) 



where U and 7 are functions of p and z. It is well known 
that these coordinates do not behave like typical cylin- 
drical coordinates. Einstein's field equations in vacuum 
with zero cosmological constant give the following linear 
partial differential equation for U, 



d 2 u 

dp 2 



18U 
p dp 



d 2 U 
dz 2 



0. 



(2) 




Equation ([2]) is Laplace's equation in a Euclidean 3-space 
in cylindrical polar coordinates. The general solution 
to ([2]) can be written out and the associated Einstein 
equations 



?l = 2 p—— (3) 

dz dp dz 



can be considered solved. Since in the weak field we have 
gu ~ —(1 + 2£7), it is tempting to consider a solution- 
generating procedure wherein one takes a known Newto- 
nian potential U (in an unphysical Euclidean 3-space), 
and solves for 7. We then have an exact solution to 
Einstein's equations. Unfortunately, the resultant Weyl 
solution bears little similarity to the Newtonian solution 
and the physical meaning of most Weyl solutions so pro- 
duced remain unclear. The culprit is the non-linearity of 
Einstein's equations that enters via ([3]). 

The Curzon Q - Chazy Q solution (CC hereafter) is 
one of the simplest special cases of the Weyl metric (JTJ) £3 
The potential is taken to be the Newtonian potential of 
a point mass (m) at the center of a fictitious Euclidean 
3-space, p — z = 0, 



U = - 



m > 0. 



(4) 
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1 For a review see 
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With fl4|) it follows from © that 
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2 2 



2(p 2 + z 2 ) 5 



(5) 



The resultant metric components are well defined except 
at (p, z) — (0, 0). The circumference of a trajectory con- 
stant t, p and z is equal to 27rp , which be- 
haves like Euclidean cylindrical coordinates for p/m or 
z/m 3> 1, but the circumference diverges in the plane 



z = as p goes to zero. Further, the meaning of m in 
the CC solution is no longer obvious. This is discussed 
below. 

Given Q and ([5]) it follows that there are two non- 
vanishing independent polynomial invariants of order p = 
2. For convenience, define 



2 2,2 

r = p + z . 



(6) 



The invariants can be taken to b^l 



wlr = 2 m ( e r 



and 



w2r = 3 m e 



(3 r 6 — 6 mr 5 + 3 m 2 r 4 + 3 m 2 p 2 r 2 — 3 p 2 m?r + p 2 m 4 ) 
r 12 (e 2 ?) 4 



(m — r) (2 r 6 — 4 mr 5 + 2 m 2 r 4 + 3 m?r 2 p 2 — 3 p 2 m 3 r + p 2 m 4 ^ 



(7) 



(8) 



Clearly, the explicit representation of (O and (j8]) is a mat- 
ter of some choice, but see below. In any event, at first 
glance, it would appear that r = is singular. However, 
along p — 0, we note that 



w2r 
~6~ 



and, in particular, 



wlr 



3/2 



m 3 (z — my 



p=0 



lim v 



0. 



(9) 



(10) 



Now, whereas \ has a local minimum (— 0) at z = ±m, 
and a local maximum at z = ±(1 ± l/^/3)m (see be- 
low), it is clear that there is no scalar polynomial sin- 
gularity along p = 0. Yet, more generally, except per- 
haps for selected trajectories, wlr and wlr both di- 
verge at p = z — 0. The directional divergence of 
K was, as far as we know, first noticed by Gautreau 
and Anderson [8j . This observation generated much fur- 
ther consideration. In terms of "polar" coordinates (r, 9) 
(p = rsin(8),z = rcos(6>), see below), Stachel @ showed 
that the area of surfaces of constant t and r decreases 
with decreasing r up to a minimum (which can be shown 



to be r ~ 0.5389m) and then diverges as r — > 0. These 
surfaces can be shown to be topologically spherical, a 
point that is important for the construction of the Hawk- 
ing mass dicussed below. Cooperstock and Junevicus [Tpj 
showed that even trajectories of the simple form z = Cp n , 
where C and n are positive constants, give K a rich struc- 
ture. Further analysis of the geodesies followed in order 
to explore the singularity and the global structure of the 
metric [ll| , [l2j , [H| , Q3 13 The general consensus is that 
the singularity has a "ring-like" , rather than "point-like" , 
structure. This is not so simple as it first sounds. The 
"ring" has finite radius but infinite circumference [l2| . 
Rather remarkably late was the computation of wlr in 
15], a work which gave visual information on the CC 
metric based on the principal null directions. This pro- 
cedure gives much less information than the visualization 
procedure considered here as we explain below. More re- 
cently, Taylor [la ] has suggested a technique for unravel- 
ling directional singularities^ 



2 Much of the older literature refers only to the Kretschmann scalar 
(K = 'R. a i3fi,v'R- a P ,n/ 1 where 1Z is the Riemann tensor) which, in 
the modern literature, is not considered a fundamental invari- 
ant. See 0] and references therein. In vacua, K = 8u>lr and so 
nothing is lost by considering K as opposed to wlr. However, 
as explained below, w2r must also be considered. Note that the 
notation wlr and w2r refers to the real part of the first and sec- 
ond Weyl invariants. Their signs are a matter of convention but 
we adhere to the conventions in [ij. 



3 It is remarkable that the simple coordinates used by Stachel 
(which we refer to as Stachel, rather than "polar" coordinates) 
have not been exploited further in the study of the CC metric. 
This is examined in Appendix A where the relevant results of 
Cooperstock and Junevicus, and Scott and Szekeres are general- 
ized. 

4 Taylor's technique is not applicable to the CC solution due to a 
critical point along p = as z — > 0. This is explained below. 
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III. GRADIENT FIELDS 

As in [l[ we define the gradient field^fl 



C = ±V a I„ = ±g af) 



din 



where now 



I\ = wlr, I2 = w2r. 



(11) 



(12) 



A. Stachel Coordinates 

The explicit forms for the gradient fields are given in 
Appendix B in Stachel coordinates. We prefer to draw 
the flows explicitly in Weyl coordinates as given below, 
but there is no loss of information in doing this since the 
inverse transformations, r = \/ p 2 + z 2 , 9 = arctan(p/z), 
are so simple^ A convention we use, based on Stachel 
coordinates, is that we say that the positive gradient flow 
is "ingoing" for k r < and "outgoing" for k r > 0. 



affect the shape of the resulting flow lines. In these plots, 
the color of the sampled flowlines is used simply to high- 
light the different regions. 




I=+3 

undefined 



B. Weyl Coordinates 

Figures [T]and [5] show the positive gradient fields for the 
CC metric in the first quadrant only, since the metric is 
axially symmetric, as well as symmetric about the equa- 
torial plane [z — > —z). The first field, ki, for the CC met- 
ric, has 7 critical points (p, z): (0, 0) (along p = Q,z 0), 
(0, ±m) (stable nodes of index +1 and isotropic critical 
points [l[ with wlr = 0), (0, ±m(l — l/\/3)) (hyperbolic 
saddle points of index —1), and (0, ±m(l + 1/V3)) (hy- 
perbolic saddle points of index —1). Note that the gra- 
dient field is undefined at (0,0) (along z = 0, p 0). 
The overall index for any hypersurface of constant t and 
4> is +1 as expectedd On the other hand, the second 
field, ki has 8 critical points: the same 7 points as fci, 
except that now (0, ±m) is an isotropic critical point with 
w2r = k L 2 = index = 0. In addition we find the critical 
points (?« 1.11m, 0) (stable node of index +1). Again we 
find that the overall index for any hypersurface of con- 
stant t and 4> is +1. Note that the field was normalized 
for visual representation, since it increases exponentially 
as we approach the origin. This normalization does not 



5 Since we are following the conventions in Q , there is a (physically 
irrelevant) reversal in the direction of the gradient flow between 
wlr and w2r. We refer to the choice " + " as a positive gradient 
flow and the choice " — " as a negative gradient flow. We find 
occasion to use both. 

6 That is, to view the flow in Stachel coordinates simply think of 
r as a circle centered on an origin at p = z = 0, and 9 a straight 
line through the origin measured from along the vertical to it/ 2 
in the equatorial plane. 

7 Of course the indices must be calculated in the full "plane" , for 
both positive and negative z and for <j> = and <j> = n. 



FIG. 1. The positive gradient field fci of the Weyl invariant 
wlr for the CC metric, presented in Weyl coordinates. The 
field is normalized for visual representation, and the flowlines 
are colored to categorize them into distinct groups according 
to their global behavior (blue, red and green). Critical points 
are represented by black circles, and critical directions of the 
fields are brown. The character of the critical points along 
with their associated indices are shown and discussed in the 
text. 
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FIG. 2. As in Figure[T] but for k 2 . 

As we can see in Figure [TJ all of the flowlines of k\ ter- 
minate at the singularity (0, 0) (approaching it along the 
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equatorial plane [z = 0,p —> 0)). They can be catego- 
rized into three distinct regions according to where they 
originate, separated by two critical directions drawn in 
brown: 

1. Blue region: The Bowlines originate from infinity. 

2. Red region: The flowlines originate from the crit- 
ical point (0,m), and the critical direction of the 
gradient field separating the blue and red regions 
originates from the critical point (0, m (1 + 1/ \/3)) • 

3. Green region: The flowlines originate form the crit- 
ical point (0, 0) along the z-axis. The critical direc- 
tion separating the red and green regions originates 
from the critical point (0, m (1 — 1 / v3)) . 

The picture for the flowlines of is somewhat dif- 
ferent, and is shown in Figure [2j However, we still have 
three distinct regions separated by two critical directions: 

1. Blue region: The flowlines originate from the criti- 
cal point ~ (1.1m, 0), and terminate at infinity. 

2. Red region: The flowlines originate from the criti- 
cal point ~ (1.1m, 0), and terminate at (0,0), ap- 
proaching it along the equatorial plane. The critical 
direction of the gradient field separating the blue 
and red regions originates from the critical point 
~ (1.1m, 0), and terminates at (0, m (1 + l/v3)j- 

3. Green region: The flowlines originate form the crit- 
ical point (0, 0) along the z-axis, and terminate at 
that point as well approaching it along the equato- 
rial plane. The critical direction separating the red 
and green regions originates from the critical point 
(0,m (1 — 1/V3)), and terminates at (0,0) . 



C. Scott - Szekeres Unfolding 

The unfolding of (p, z) = (0,0), given by Scott and 
Szekeres in [f| and [13j . is reproduced in Appendix C. 
Spacelike infinity (r — > oo) is mapped onto < X < tt 
with Y = tt/2 and X — tt with < Y < tt/2, the z axis 
maps onto — tt/2 < Y < tt/2 with X = and the p axis 
onto 7r/2 < X < tt with Y = 0. Now r = corresponds 
to < X < tt/2 with Y = -tt/2 and X = vr/2 with 
—tt/2 < Y < 0. The singularity of the CC metric is 
represented only by X = tt/2 with Y = 0. 

The gradient fields of the CC metric, with this un- 
folding, are shown in Figures [3] and |U The unfold- 
ing of the Weyl coordinate point (0, 0) is now evident: 
All of the flowlines for k\ terminate at the singular- 
ity (X, Y) — (7r/2,0). The flowlines of the green re- 
gion now originate from a different point, specifically 
(X,Y) = (tt/2, -tt/2). 




FIG. 3. As in Figure [T] but with the Scott - Szekeres 
unfolding. All flow lines terminate at the singularity at 
(X, Y) = (tt/2,0). All green flow lines originate at (X, Y) = 
[tt/2, — tt/2) which is the critical point (p, z) = (0,0) along 
p = 0,z -> 0. 




FIG. 4. As in Figure[2]but with the Scott - Szekeres unfolding. 



D. A New Unfolding 

Whereas the unfolding of (p, z) = (0, 0) given by Scott 
and Szekeres accomplishes the task, the rather compli- 
cated procedure also modifies the entire spacetime repre- 
sentation. Here we seek a new unfolding of (p, z) = (0, 0) 
which does not modify the spacetime in the large. It 
turns out that we need only modify p. On reviewing Ap- 
pendix A we see that the most important term to consider 
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is the exponent of sm(9)/x. We can write 
sin(0) 



pra 



(13) 



Now either this exponent diverges or it does not. If it 
diverges we can compactify this divergence with a tanh 
function. If it does not diverge we can set the term to 
zero by multiplying by p. Finally, let us require that the 
new p and old p approach each other for sufficiently large 
r. We arrive at the unfolding 



m to \m / V to 



(14) 



The singularity is at (p, 0) = (m, 0) and the critical point 
is at (p, 0) = (0,0). No flow lines cross the "edge" < 
p < m where r = 0. Indeed, whereas the trajectories 9 = 
and 9 — tt reach p = 0, all other trajectories of constant 
9 reach p = to. This is shown in figure [SJ The gradient 
fields of the CC metric are shown in Figures [6] and [7] 
with this new unfolding. It is very important to realize 
that this unfolding cannot correct all misrepresentations 
created by the Weyl coordinates. In particular, if the 
Weyl coordinates do not cover the "edge", neither does 
the unfolding. 
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FIG. 6. As in Figure [T] but with the new unfolding. 
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FIG. 7. As in Figure [2] but with the new unfolding 



FIG. 5. The z/m — p/m quarter plane. These are not oblate 
spheroidal coordinates. The "edge" < p < m, where r — 0, 
is not part of the spacetime. 



E. "mass" 

Whereas the constant m has entered the CC solution 
via the Newtonian potential (j4|), the meaning of to within 
the CC solution is no longer obvious. This is explored 
in Appendix D where we show that to is certainly the 
"mass" at spatial infinity. However, away from spatial 



infinity, looking at quasi local and local constructions, 
we find that the Hawking mass AAh provides no useful 
information for the CC solution. Rather, it is the classical 
effective gravitational mass M. = lZ g ® ^g g ^ 2 /2 [l| that 
provides useful information on the CC solution. Whereas 
this M. rapidly converges to to with increasing r, near 
r = 0, A4 shows considerable structure. Most interest- 
ing is the fact that A4 = at the (naked) singularity, 
reminiscent of spherically symmetric naked singularities 
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IV. THE NEWTONIAN ANALOGUE 

We finally get to the central consideration of this pa- 
per. As explained previously following for example 
Ellis [2l|, we construct the Newtonian tidal tensor 



1 



,a,b ~ ^VabO 



E ab = <?> 

and define the associated invariant 



E ab E ab = $ a!6 $> a > fe -i(n $) 2 . 

We now construct the gradient field 
l\ = ±V c {E ab E ah ). 



(15) 



(16) 



(17) 



We say that l\ is a Newtonian analogue (in no way any 
limit) of ki if their associated phase portraits are "anal- 
ogous" , a designation which is explained in detail be- 
low. This analogy is strict since, as explained in [![, 
kf = ± V° {Ep-yEP~<) where Ep y is the usual electric 
component of the Weyl tensor; that is, the general - rel- 
ativistic tidal tensor in the Ricci - flat case. The gener- 
alization of (|17l) for comparison with k% is 



1^ = ±V d (E b a EZE^). 



(18) 



The Newtonian potential for a infinitely thin ring in 
vacuum with radius a and mass m is 



ring 



271 



dO 



yjp 2 + z 2 + 1 - 2pcos6» 



(19) 



where x = x/a Vi. Note that m only affects the inten- 
sity of the gradient field, but does not change the shape 
of the flowlines or the normalized field when the coor- 
dinates are parameterized by the radius. We therefore 
take m = 1. Further, since there is no scale set for z, 
we take y/2z for visualization. Along the z-axis, for both 
1 1 and I2, defined by p7|) and (Tig)) respectively, we find 
5 critical points: s/2z/a = 0, ±1, and ±\/3 (~ ±1.7). 
This is to be compared to : z/m — 0, ±1, ±(1 — 1/V3), 
and ±(1 + 1/V3) (~ ±1.6) for the CC metric. The New- 
tonian gradient fields are shown in Figures and [9l In 
both cases the overall index is +1. 

Comparing Figures [6] and [8] we see that the CC metric 
contains an additional flow (the green flow) not present 
for a Newtonian ring. Otherwise, the gradient flows are 
remarkably similar. Comparing Figures [7] and |H] we see 
that as well as the additional flow in the CC metric, there 
is also an additional critical point in the equatorial plane 
not present for a Newtonian ring. Again, otherwise, the 
gradient flows are remarkably similar. 

The green region in the CC solution, which is absent 
in our Newtonian analogue, indicates that the "source" 
of the CC metric is not just a Newtonian ring. 




FIG. 8. Gradient field h, defined by (fT7|l . for a Newtonian 
ring in vacuum. Compare Figure [6] 



42 2 




TBI), for a Newtonian 



FIG. 9. Gradient field l 2 , defined by 
ring in vacuum. Compare Figure [3 



V. DISCUSSION AND CONCLUSION 

We have demonstrated a new tool designed to visualize 
spacetime curvature based on the construction of gradi- 
ent flows of invariants. The emphasis here has been on 
the construction of strict Newtonian analogues wherein 
the invariants are the relativistic and Newtonian tidal 
invariants. The case we have considered, the Curzon - 
Chazy solution, is by no means an easy case to start 
out with. Despite this, we have shown that the gradi- 
ent flows for the tidal invariant in the Curzon - Chazy 
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solution, and that for an infinitely thin Newtonian ring, 
are remarkably similar in detail. Only "interior to" and 
"closely above" the Curzon - Chazy "ring" do the flows 
differ. This region is absent in the Newtonian case. For 
completeness we have included all non - differential in- 
variants of the Curzon - Chazy solution and so we have 
also included a study of the second Weyl invariant and its 
Newtonian counterpart for an infinitely thin ring. Again 
we find that the gradient flows are remarkably similar 
in detail except for the same region close to the Cur- 
zon - Chazy ring and the addition of an external, but 
"close" , critical "ring" of the flow in the equatorial plane 
not present for the Newtonian ring. Whereas one could, 



perhaps, refine the Newtonian counterpart and, perhaps, 
find sharper agreement, the arguments presented here, 
we believe, go a long way to clarify the notion that the 
"source" of the Curzon - Chazy solution is "ring - like" 
and that the construction of strict Newtonian analogues, 
correct "in the large" , is possible. 
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Appendix A: Stachel Coordinates 



Under the transformations 

p = rs'm(6), z = rcos(0) (Al) 

the CC metric takes the form 

ds 2 = -e-^dt 2 + {e- m2 7i " )2 (dr 2 + r 2 d6 2 ) + r 2 sin((9) 2 d</> 2 ). (A2) 
Either from (|A2p , or from CC metric and (|A1|) . we obtain 

wl = (wlr)m = 2e ^ (A3) 



and 



„ , „ \ « „ vi {x~ I)(2t2 + sin(6») 2 t3) 
w2 = (w2r)m 6 = -3e 3tl - T — - (A4) 
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where 

x = r/m, (A5) 
sin(9) 2 -2x 

tl = , (A6) 

x z 

t2 = x 2 {x-l) 2 , (A7) 

and 

t3 = 3x 2 - 3x+ 1. (A8) 

Now let us take 

sin(0) = a(x), (A9) 
where < a < 1, in order to define trajectories 9(x). We are interested in limits as x — > 0. First let us assume a £ C 2 . 
i) a(0) = a D ± 
We find 

a 2 fl 2 

wl ~ 2e 2 ^-£r -> oo, (A10) 
a; iu 

and 

a 2 „2 

w2 ~ -3e 3 ^ -> -oo. (All) 

nj o(0) = 
We find 



and 



wl ~ -> 0, (A12) 

e» a; 8 



w2 I^H ^ °- (A13) 



Historically, a € C° have played a role. Following 10], in the notation of [l2j], consider the trajectories 

-=b(bH-T, (b,n>0), (A14) 
m 2 m 



that is, 



xcos{6) = b(-)i(xsm{9)) n . (A15) 



Now if a(0) = 0, we can use the small angle formula to obtain 



o=^jp S = 6(i)i (A16) 



Of particular interest is the case n = 2/3 for which 
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For the convergence of wl and w2 we must have 

a 2 - 2x 



< 0, (A18) 

that is, b > 1 for the case (|A17[) . This is the correction in [l2[ to an error in Further details concerning the 

very particular choice (|A14[) can be found in 11211 . Here we simply note that a £ C 1 requires n < 1/2 and a G C 2 
requires n < 1/3 for this particular choice. In [1J], Scott and Szekeres go on to argue (on page 562) that they found a 
trajectory along which the Kretschmann scalar goes to a finite non - zero constant as r — ¥ 0. We have examined this 
claim in detail and find the claim to be false. Along the suggested trajectory we find that the Kretschmann scalar 
goes to zerojfl Our understanding is that only two limits are possible within known coordinates: zero and ± infinity. 



Now for the metric (|A2[) . £ Q = 5" is a Killing vector. As a result, for all geodesies with tangents u a we have a 
constant of the motion £, a u a = —7. As a result, for all geodesies we have 

dt 2 . 

— = ± 7 e" (A19) 

where A is an affine parameter. This shows the inadequacy of the coordinate t as x — > 0. 

Appendix B: Gradient Fields in Stachel Coordinates 

The contravariant gradient flow associated with the invariant wlr is given by 

= ±2e3tl sin W2 (4sin(^3 + M) + 6t5 

and 

^ =T4e «iM^(g)(yfl+«i) (B2) 

where 

i4 = x(36x 3 - 63x 2 + 34a;-4), (B3) 
t5 = x 3 (3x 2 -6x + 2)(x-l), (B4) 

and 

t6 = x 2 (9x 2 -15a; + 7). (B5) 

The contravariant gradient flow associated with the invariant w2r is given by 

7 m 6(x - l)t5 + sinC^^esin^) 2 ^ - l)t3 + xt7) 
m «2 = T3e ^ (B6) 

where 

<7 = 45x 4 - 126x 3 + 124x 2 - 50a; + 6, (B7) 

and 

m 8 k e = ±6c 4ti Sin(g) cos(e)(x - 1)(3 sin(g) 2 ^3 + t6) 
2 x ls 

Note that the upper sign corresponds to a negative gradient flow in all cases above. 



8 This is rather unfortunate in the sense that had the claim been necessary for a complete understanding of the singularity, 

correct, a further unfolding of the singularity would be absolutely 
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Appendix C: Scott - Szekeres Unfolding 

The unfolding of (p, z) = (0, 0) given by Scott and Szekeres in [l2| and [13| . obtained by trial and error, is 

X = arctan((p/m) e m / z ) + arctan((p/m) el-^/p)"*) (CI) 

and 

F = arctan ^3- - ((r/m)8 + : + i (p/m)2 (r/m) - 4) i/4 J ( C2 ) 
in Weyl coordinates. In Stachel coordinates we find 

X = arctan^sin^e 1 /*"^)) + arctan^sin^e^^ 8 ^^ 273 ) (C3) 

and 

Y = arctan 3x cos (0) - v , In . C4 

V (a; 8 + 1 + sm(6») 2 /3a; 2 ) 1 / 4 y 

Appendix D: "mass" 

We are interested to see how the constant m in the CC solution is related to "mass". First, let us rewrite (|A2[) via 
Taylor series about 1/r = with explicit terms to order 1/r. We have 

ds 2 = -(1 - ^)dt 2 + (1 + ^)(dr 2 + r 2 dn 2 ) (Dl) 

where dQ 2 is the metric of a unit 2-sphere (d9 2 + sin 2 (9) dip 2 ). To this order in r we find that the Einstein tensor for 
(ID1|) vanishes. Due to the spherical symmetry of (|D1[) we can consider the mass defined by M. = lZ e ® ^g g ^ 2 /2 [l|. 
To order 1/r we find 

3m 2 /1\ , s 

M=m-—^-j. (D2) 

As a result, at spatial infinity, m is the "mass" . We need not consider the ADM mass, the Komar integrals nor the 
Bondi -Sachs mass [13] • Rather, we are interested in exploring the meaning of m away from spatial infinity. We are, 
therefore, interested in local and quasi-local quantities. 

Let us look at the Hawking mass [18| , which can be defined by 

M "=(#)" 2 ( 1 - tJ/<^) < D3 > 

where S is a spacelike topological two-sphere, A(S) is the associated area and here p and (x are the Newman - Penrose 
spin coefficients. As mentioned above, Stachel [9( showed that for S defined by surfaces of constant t and r in (|A2[) . 
A decreases with decreasing r up to a minimum (which we find to be r ~ 0.5389m) and then diverges as r — > 0. To 
ensure that S is a topological two - sphere we use the standard Gauss - Bonnet theorem 



Ky/gdx a dx b = 2n X (S) (D4) 



where x a are the coordinates on S, g is the determinant of the metric on S, 1Z is the Ricci scalar on 5, and x(^) i s 
the Euler characteristic of S. Again taking S to be defined by surfaces of constant t and r in (IA2I) we find 



X (S) = 2, (D5) 

and so we are indeed considering topological two - spheres. Constructing a complex null tetrad in the usual way, we 
find that for El 
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and so we find 
M 



^ = -i/32 V2 e yJ r 0§==d6 f-a + (* e-(?-y+^<»fy)y (2 y - 2 - y 2 + (cos W) 2 y 2 ) 2 <w) p- 1 

TO V Jo Ve( sm ( e )) » V Jo v y / 



(D7) 



where y = m/r = 1/x. Expanding about y = we find that 

.M/f 37r 1/1 tt\ r . m ,m 



4 2 \ 2 4/rn r 



and so, for large r, A4h is not a reasonable measure of "mass" for the CC solution. Examining (|D7[) in more detail we 
find that A4h = for r/m ~ 8.86, .M// reaches a maximum Ain/fn ~ 1.25 at r/m ~ 2.2 and — > — oo as r — » 0. 
We have to conclude that A4h is nowhere a measure of "mass" for the CC solution. In this regard, one is reminded 
of the result of Hansevi [l!| who showed that the Hawking mass can be negative even for convex two - surfaces in 
Minkowski spacetime. 

Returning to M = lZ g ® ^g gg 3 ^ 2 /2, whereas this quantity is usually restricted to strict spherical symmetry, it is 
well defined for (|A2[) . Indeed, a straightforward calculation gives 

M _a f tl 



= e 



m \ 2 

where tl is given by (|A6|) . Continuing with (|A9I) we find 



1 + TT ), (D9) 



for a(0) = a a ^ and 



lim ( — ) 0, (D10) 



lim f — ] ->■ -oo, (Dll) 
x->0 \ m 



for a(0) = a Q = 0. In strict spherical symmetry, it is known that naked singularities are massless (Ai = 0) [20]. The 
result (|D10[) suggests that this might hold away from spherical symmetry. It should be noted that — converges very 
rapidly to 1 for all 8 with increasing x. We conclude that At, if not the "mass", at least summarizes some important 
properties of the CC solution. 

I 



